Extension of the full-folding optical model for nucleon-nucleus 
scattering with applications up to 1.5 GeV 

H. F. Arellano 

Departamento de Fisica, Facultad de Ciencias Fisicas y Matemdticas 
Universidad de Chile, Casilla 487-3, Santiago, Chile 

H. V. von Geramb 

^ \ Theoretische Kernphysik, Universitat Hamburg 

Lumper Chaussee 149, D-22761, Hamburg, Germany 

<' 

in 

Abstract 



O 



o 

=5 



X 



(Dated: February 8, 2008) 



The nonrelativistic full-folding optical model approach for nucleon-nucleus scattering is extended 
into the relativistic regime. In doing so, kinematical issues involving the off-shell Lorentz boost 

o 

of the colliding particles between the two nucleons and the projectile-nucleus center-of-mass ref- 
erence frames have been taken into account. The two-body effective interaction is obtained in 



the framework of the nuclear matter ^-matrix using nucleon-nucleon optical model potentials that 
fully account for the inelasticities and isobar resonances in the continuum at nucleon energies up 
to 3 GeV. Diverse nucleon-nucleon potential models were constructed by supplementing the basic 
Paris, Nijmegen, Argonne or Gel'fand-Levitan-Marchenko inversion potentials with complex sep- 
arable terms. In each case the additional separable terms ensured that the combination led to 
NN scattering phase shifts in excellent agreement with experimental values. With each phase shift 
fitting potential nuclear matter ^-matrices have been formed and with each of those relativistic 
full-folding optical potentials for nucleon-nucleus elastic scattering determined. Application to such 
scattering for projectile energies to 1.5 GeV have been made. Good and systematic agreement is 
obtained between the calculated and measured observables, both differential and integrated quan- 
tities, over the whole energy range of our study. A moderate sensitivity to off-shell effects in the 
differential scattering observables also is observed. 

PACS numbers: 03.65.Nk, 25.40.Cm, 24.10.-i 
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I. INTRODUCTION 



Elastic nucleon-nucleus (NA) scattering is known now as an excellent means of testing 
nuclear structure Q. Results of the nonrelativistic theory for such scattering confidently 
can be used as predictive for information required in applied nuclear technology, where 
large amounts of nuclear reaction data -including fission cross sections at intermediate 
energies- are required for various challenging applications such as accelerator transmuta- 
tion of waste, particularly the elimination of long-lived radioactive wastes with a spallation 
source, accelerator-based conversion to destroy weapon grade plutonium, accelerator-driven 
energy production to derive fission energy from thorium with concurrent destruction of the 
long-lived waste and without the production of nuclear weapon material, and accelerator 
production of tritium ||. There is also a great need for such information to be the base in 
analyzes of patient radiation therapy and protection. 

With basic science there is the intellectual challenge to go beyond the physics of a single 
hadron and understand essential aspects of nuclear physics from first principles such as QCD 
. It is generally agreed that the QCD Lagrangian involves nonlinear dynamics. Such makes 
it very difficult then to understand nuclear physics fully from first principles, and so most 
of nuclear physics phenomena are interpreted in terms of appropriate effective degrees of 
freedom. One such view of natural phenomena in terms of energy scales (Q), divides the 
nuclear-hadronic scale into the nuclear structure region Q ~ 1 — 10 MeV and the nucleon 
and nucleon-nucleon (NN) region, with structure and substructure scales Q ~ 0.3 — 1 GeV. 
This large separation between the hadronic energy scale and the nuclear binding scale poses 
stringent difficulties to apply nonlinear QCD directly to understand the physics of nuclei. 
However, quantitative calculations based on effective quantum field theory (EQF) techniques 
that arise from chiral symmetry provide an alternative approach. At present, this method 
is being extended to address few- and many-nucleon interactions. When combined with 
first principle calculations of the low energy constants from QCD, these EQF may provide 
a consistent qualitative understanding of properties of nuclei and of low to medium energy 
nuclear scattering. 

Besides understanding the structure of nuclei from a QCD point of view, it is of interest 
to understand the dynamical behavior of nucleons in the presence of nuclear matter. The 
relevance of modified NN scattering amplitudes in the form of ^-matrices in the nuclear 
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medium with mean fields and Pauli blocking is well known. These NN amplitudes and 
alternative reductions -in the form of t-matrices- have been used with qualitative success 
in the specification of nuclear densities in stable nuclei and the description of NA scattering 
for projectile energies (T^) below 1 GeV. For above this limit we expect significant 
dynamical changes due to dibaryonic fusion with subsequent fission in the short range region 
of NN subsystems J§ . 

Hitherto quite independently, several groups successfully described intermediate energy 
NA scattering using two nucleon t- or g-matrices as driving effective interactions ||. Among 
them we distinguish two main philosophies under which a description of NA collisions is made 
and which we specify as the nonrelativistic Schrodinger approach and as the relativistic 
Dirac approach Jl], |5fl. Common to both is the explicit use, and accurate treatment, of 
an interacting NN pair in the realm of other nucleons in a nucleus, and the need for an 
effective NN scattering amplitude known and defined on- and off-shell. For the more recent 
calculations of nonrelativistic NA optical potentials in momentum space || [j], § , a consistent 
treatment of the fully off-shell NN t- or g-matrices have been used. 



As nearly all available NN potentials ||, 1C, 11, 12, 13] have been fitted to NN scattering 
data up to 350 MeV, a nonrelativistic treatment of the full-folding optical model suffices 
most NA scattering applications consistent with the input NN model. Recently, however, 
inversion potentials developed to fit elastic NN phase-shifts for T Lab < 1.2 GeV have been 
used within the full- folding approach for projectile energies as high as 500 MeV . De- 
spite the fact that these were nonrelativistic applications, they allowed a better description 
of data relative to what was found previously using traditional NN potential models. The 
lack of relativistic kinematics as well as pion production and A excitations within the NN 
pair did not manifest itself as a dramatic limitation of the model, even though the theoretical 
confidence level was reached. It is the primary aim of this work to remove these shortcom- 
ings and formulate a momentum space full-folding model using minimum- relativity. With 
the appropriate modifications to the existing analysis programs we expect to extend the 
confidence level into the GeV region and to obtain high quality full-folding optical model 
results at energies as high as 1.5 GeV. 

Extensive studies of NA scattering in the context of relativistic Dirac models have been 
made j5|, [T5], [T^, [T7|, While these models are closer to a fundamental formulation with 
the inclusion of relativistic kinematics and dynamics within a single framework, to date they 
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all have used fitted NN interactions that do not describe the phase shift data well. This is 
a serious drawback as it is well established that any folding model requires a high quality 
NN interaction as input. 

The description of NA scattering requires only modest nuclear matter densities and thus 
the need of NN ^-matrices is less demanding. Furthermore only for projectile energies 
below 500 MeV are nuclear medium effects predominant with specific treatment of Pauli 
blocking and self consistent mean field effects being crucial. For medium and higher NN 
energies, meson production and intrinsic hadron excitations within the interacting pair are 
important. Consequently the t-matrix no longer is unitary in the elastic channel. However as 
the low energy t- and ^-matrices are well defined within potential scattering theory, we seek 
a continuation of NN potential models with an NN optical model potential (NNOMP). We 
have devised and generated NNOMP for 0.3 < < 3 GeV and applied them to NA elastic 
scattering. The calculations are based upon a relativistically corrected full-folding optical 
model in momentum space that is an extension of a nonrelativistic predecessor || . Since the 
Lorentz contraction scales as the ratio of the projectile energy to its mass, it requires nucleon 
projectile energies above 400 MeV to have sizeable contraction effects. It is safe to include 
relativistic kinematics in folding calculations when the projectile energy surpasses 300 MeV. 
Relativistic kinematics is widely used in pion- nucleus scattering |B|, [20|, |2T], ^ |23| . An 
excellent discussion on some aspects about relativistic kinematics in NA scattering is given 
in the review by Ray, Hoffmann and Coker ||. 

One of the important advantages of finding a nonrelativistic NA optical model potential 
is its well defined structure in terms of interacting nucleon pairs. This framework has been 
remarkably successful in the study of low energy scattering with its link to the nuclear 
shell model of single particle bound states. The combination of target correlations and high 
quality NN interactions has provided a good first order description of the NA dynamics. 
This nonrelativistic theoretical approach has been applied with success for many energies 
and targets and now it is timely to incorporate relativistic corrections so that calculations 
can be made with energies as high as 1.5 GeV. 

A brief summary of the current full-folding optical model is given in Section II and in 
which the points where relativistic kinematic corrections are to be made are stressed. The 
current experimental situation of NN phase shift analysis up to 3 GeV is discussed in Sec- 
tion III. We describe also how the NNOMP are determined by the NN phase shift data. 
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In Section IV we set the framework of the g-matrix calculations and specify minimal rela- 
tivity considerations. In Section V some details of the full-folding calculations and various 
applications are presented. We discuss the role of the Fermi motion in the NN effective inter- 
action and analyze the sensitivity of NA scattering observables upon the use of alternative 
approaches for the relativistic kinematics. We also examine total cross sections for nucleon 
elastic scattering and differential observables at beam energies from the hundreds of MeV 
up into the GeV regime. Section VI contains a summary and the main conclusions of our 
study. Finally, we have included two Appendices. In Appendix A the relativistic kinematics 
transformations in the context of the full-folding approach are outlined as extracted from the 



articles by Aaron, Amado and Young [21] and by Giebink [23]. In Appendix B we present 
the algorithm we have used to determine the NNOMP from data. 



II. FULL-FOLDING FRAMEWORK 



In the nonrelativistic theory of the optical model potential, the coupling between the 
projectile and the target nucleus in the elastic channel is given by the convolution between a 
two-body effective interaction and the target ground-state mixed density. In the projectile- 
nucleus cm. frame, the collision of a projectile of kinetic energy E is described by the optical 
potential U(E) which in a momentum space representation is expressed as || 



U(k>, k; E) = Y, J J dp'dp<^pO<ky|T(n a )|kpWa(p) (i) 

where (fi a represents target ground state single-particle wave functions of energy e a , and 
a < tF restricts the sum to all levels up to the Fermi surface €f- The two-body T matrix is 
evaluated at starting energies Q a = m p + E + m t + e a consisting of the mass of the projectile 
m p , its kinetic energy E, the mass of the target nucleon m t and its binding energy e a . The 
subscript A+l indicates matrix elements in the projectile-nucleus cm. frame and recoil 
effects have been neglected for simplicity. In the most recent full-folding calculations of the 
optical potential for NA scattering, where explicit medium effects are incorporated in the 
two-body effective interaction, the T-matrix is represented by an infinite nuclear matter 
g-matrix ||. In the absence of these medium modifications, the T-matrix has usually been 
approximated by the two-body scattering matrix associated with the collisions of two free 
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particles |], [7|. None of these full-folding approaches incorporate the necessary relativistic 
kinematics needed for high energy processes. 

The above expression for the optical potential requires the T matrix in the projectile- 
nucleus cm. frame. However, most practical two-body potential models are designed to 
account for the scattering data in the two-body cm., where a one-body wave equation 
(a Schrodinger kind equation) is used to construct the realistic two-body bare potential 
fit to elastic scattering and ground state data. The practical problem that emerges then 
is how to make use of this description to extract the needed effective interaction in the 
projectile-nucleus cm. with an adequate account of relativistic effects. This has been 
a long standing problem in nuclear research and various approaches have been proposed 
and discussed elsewhere [19], [21], ^2], p3| . In the procedure followed here we retain the 
dynamical structure of the optical potential as expressed by Eq. ([!]) and identify, via a 
Lorentz boost, the corresponding kinematical variables involved in the two-body collision. 
Thus, the transformation of the T matrix from the two-body (2B) to the projectile-nucleus 
(A+l) cm. frame can be done by considering three separate aspects. First, as Lorentz 
invariance of the flux is required, an overall normalization factor -usually referred as the 
M0ller flux factor- mediates between the scattering amplitudes in the two frames. Second, 
the kinematics in the projectile-nucleus cm. frame needs to be transformed to the two- 
body cm. system, which is the reference frame where the bare potential model is defined. 
And third, the transformation of the scattering matrix from the two-body to the projectile- 
nucleus cm. frames involves the rotation of the spins, an effect referred as the Wigner 
rotation. This contribution has been studied in the context of the relativistic 'no-pair' 



potential for nucleon- nucleus scattering by Tjon and Wallace [L7| and was observed to yield 
rather moderate effects at nucleon energies between 200 and 500 MeV. Although there is 
no statement about the importance of Wigner rotation contribution at the higher energies 
considered here, we shall neglect them in the present work. 

Consistent with the notation introduced in Fig. [I] for the coupling between the pro- 
jectile and the target struck nucleon, we denote their respective incoming and outgoing 
four-momenta 



> 



k = (cu,k), k' = (u',k' 
p = (e,p), p' = (e',p') . (2) 
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When translational invariance is assumed, such is the case of free or particles in infinite 
nuclear matter, the two-body interaction is characterized by the starting energy Q and the 
total momentum of the colliding particles Q, 

(ky|T(fi)|kp> 4+1 = 77(k'p';kp) <k r ',-k r '|r Q (v^)|k r ,-k r ) 2B 8(Q' - Q) . (3) 

Here the Dirac delta-function makes explicit the total three-momentum conservation of the 
two-body collision, Q' = Q, where 



Q = k + p , Q' = k' + p', 

and -^/i represents the energy in two-body cm., 

s = n 2 - q 2 . 

The overall coefficient r\ is the M0ller flux factor, 

~u(k')e(-k')u(k r )e(-k r 



rj(kp ; k p) 



1 1/2 



(4) 



(5) 



(6) 



co(k') e(p>) w(k) e(p) 

where the energies 10 and e are on-mass-shell, i.e. u 2 = m 2 + k r 2 and e 2 = m 2 + k r 2 . 

What remains to be specified is a Lorentz transformation for the relative momenta k r and 
k r '. To this purpose we followed the approaches introduced by Aaron, Amado and Young 
(AAY) [El] and Giebink [E3J. Although, both assume a representation of the scattering 
matrix in the two-body cm., they differ in the way the Lorentz boost is devised. Details are 
given in Appendix A where we show that in both cases the relative momenta can be cast 
into the form 



k r ' 



Wk- (1- W)p , 
W'k' -(l-W')p' 



(7) 



with W and W being functions of the momenta of the colliding particles (Eqs. flA4| ) or 
( |A11|) ). At low beam energies both prescriptions meet the nonrelativistic limit W ~ W ~ 

m p / {m p + m t ), as expected. 



III. THE NN OPTICAL MODEL POTENTIAL 



Of the whole spectrum, low and medium energy NN scattering traditionally is described in 
terms of Hermitian potentials. At medium energies, production processes and inelasticities 
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become possible and several elementary systems composed of nucleons and mesons contribute 
to NN scattering. At present there is no high quality description of NN scattering above the 
inelastic threshold either in terms of QCD or in terms of nucleons and mesons. 

A high quality fit of on-shell t-matrices by means of a potential model is very desirable as 
it provides extensions of the effective interaction into the off-shell domain and into a nuclear 
medium, which are important dynamical features in few and many body calculations. Many 
examples using microscopic optical model potentials for elastic nucleon-nucleus scattering 
and bremsstrahlung reactions, have shown that it is crucial to have on-shell t-matrices in 
the best possible agreement with NN data at all energies. Concomitantly one needs high 
precision NN data against which one can specify NN interactions. To this purpose we 
have relied on a large body of experimental NN data whose parameterization in terms of 



amplitudes and phase shifts are smooth for energies to 3 GeV [f24| . This is a supposition for 
the construction of an NN potential above 300 MeV. 

There are many studies of few and many body problems in the low energy regime < 
300 MeV and the results have consequences for any model extension above threshold. We 
note in this context that significant off-shell differences in t-matrices are known to exist 
among the theoretically well motivated boson exchange models of NN scattering. It remains 
difficult to attribute with certainty any particular dynamical or kinematical feature with 
those differences. Non-locality, explicit energy dependence and features associated with 
relativistic kinematics are some possibilities. In contrast, there is the quantum inverse 
scattering approach by which on-shell t-matrices can be continued into the off-shell domain. 
A specific method is the Gel'fand-Levitan-Marchenko inversion algorithm for Sturm-Liouville 
equations. This approach to specify t-matrices off-shell is appropriate when the physical S- 
matrix is unitary and the equation of motion is of the Sturm-Liouville type. Such is valid 
without modification for NN t-matrices in the energy regime below 300 MeV, and for the 
unitary part of the S-matrix above that energy. 

In the spirit of general inverse problems, we have extended the available low energy poten- 
tial by additional complex potential terms which are determined from a perfect reproduction 
of the experimental data (here the partial wave phase shift analysis) for all energies above 
300 MeV. By that means NN optical models were generated, separately for each partial 
wave. The algorithm we have developed allows studies of complex local and/or separable 
potentials in combination with any background reference potential 0, [25|]. Here we limit 
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the reference potential to the well known real coordinate space potentials from Paris ||, 
Nijmegen Jll| (Reid93, Nijmegen-I, Nijmegen-II), Argonne [12|] (AV18), and from inversion 
p6| . To them we add channel-dependent complex separable potentials with energy depen- 
dent strengths [||. For a given input data set there is a unique NNOMP within a given 
potential class. Some of these issues are outlined in Appendix B. 

NN scattering is a long standing problem which has been reviewed often as the database 
developed. The low energy data has been analyzed by the VPI/GWU group |24| for J^ab < 
400 MeV, the Nijmegen group with the NN phase shift results PWA93 for T Lab < 350 MeV, 
and by Machleidt |L3| giving the Bonn-CD-2000. Of these, the VPI/GWU group has given 
many solutions over the years, the latest are for energies to 3 GeV ^7j. We have used 
the solutions SP00, FA00 and WI00 in our calculations and found results that differ but 
marginally. Thus hereafter in the main we refer solely to the results of calculations based 
upon SP00. 

As with the NN phase shift analysis, one boson exchange potentials have received several 
critical reviews fl3fl , including observations that there are small variations between phase 
shift analysis and potential models below the subthreshold domain Tiab < 300 MeV. A 
theoretically stable result would require many quantities, that need be specified a priori, 
to be determined by independent sources. At present that does not seem feasible and all 
current potentials rely upon fits of many of their parameters to the same data. All such 
fits, however, have been made independently of each other and are based upon differing 
theoretical specifications of the boson exchange model dynamics. 

Above 300 MeV, reaction channels open and the elastic channel S-matrix no longer is 
unitary. Only the A(1232) resonance has a low energy threshold and a relative small width 
of 120 MeV. Therefore it is the only resonance we expect to be obviously visible in the energy 
variation of the elastic scattering phase shifts. In particular one notices typical variations 
in the l D2, 3 i 7 3, and 3 Pi<2 channels. Otherwise the phase shifts to 3 GeV vary smoothly 
as functions of energy. Together with the strong spin-isospin coupling, this property infers 
optical potentials that are channel dependent in contrast to the NA case for which assumed 
central and spin-orbit potentials are partial wave independent. The plethora of reaction 
channels that open to 3 GeV, and the requirement of an NN optical potential prescription, 
mean that it is an interesting task for a microscopic model to link QCD substructures to 
NN scattering phase shift functions in analogy to that successful prescription by which NA 
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optical potentials have been determined by folding effective interactions. 

To describe this developing system for 0.3 < < 3 GeV we used Feshbach theory to 
specify the optical potential. An important feature of that theory is the projection operator 
formalism with P and Q subspaces, which divide the complete Hilbert space (P + Q) into 
the elastic scattering channel, the P space, and the inelastic and reaction channels, the 
Q space. This infers a complex and separable component in the optical potential with an 
energy dependent strength. If a very large number of intermediate states contribute, the 
effect equates to a local potential operator. 

A covariant description of NN scattering formally is given by the Bethe-Salpeter equation 

M = V + VGM , (8) 

where Ai are invariant amplitudes that are based upon all connected two particle irreducible 
diagrams. This equation serves generally as an ansatz for approximations. Of those, the 
three dimensional Blankenbecler-Sugar reduction is popular and sufficient for our purpose 
to define an NN potential |T3] . The amplitudes are now expressed with the reduced terms 
and they satisfy a three-dimensional equation 

f d 3 k M 2 AtJk)A±J-k) 

Taking matrix elements with only positive energy spinors, an equation with minimum rela- 
tivity results for the NN t- matrix, namely 

/d 3 k M 2 1 

W V(q ' k) ^i^T¥ T(k ' q) - < 10) 

Using the substitutions 
and 



nq ', q)= (g 5 v W , q) g) 3 , ( 12 ) 

a simplified form of the t-matrix is obtained. It is the familiar Lippmann-Schwinger equation 

/H 3 k M 
(^F^'^q^rk^^'^- (13) 
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The strategic importance of this result is that it defines a sensible continuation of a t- 
matrix, constrained on-shell by the experimental data and phase shift analysis, into the 
off-shell domain as required by the full-folding optical model. Thus we do not rely primarily 
on a fundamental theoretical result but rather on experimental NN data and the moderate 
sensitivity of NA scattering to alternative off-shell continuations. It is an important result 
also of this analysis that on-shell equivalent NN optical model potentials yield very similar 
NA scattering observables - irrespectively of differences in the off-shell domain and the 
constraint off-shell continuation defined with Eqs. (0-13). 



IV. IN MEDIUM EFFECTIVE INTERACTION 

A crucial step in the description of NA scattering processes has been the definition of an 
effective interaction based on the bare two-body interaction in free space. Such has been 
the philosophy of the early (local) folding models and the most recent nonrelativistic full- 
folding models 0. The effective interaction, T in Eq. ([[]), is obtained in the framework 
of Brueckner-Bethe-Goldstone for the ^-matrix. The extension of this approach to high 
energy applications requires a minimal account of relativistic corrections. Along this line we 



have followed the discussion by Brockmann and Machleidt [p8~ |, where a relativistic three- 
dimensional reduction of the Bethe-Salpeter equation is used to describe the interaction 
between nucleons in the nuclear medium. If only matrix elements between positive-energy 
spinors are taken, then the medium-modified invariant amplitude in an arbitrary frame reads 
(cf. Eqs. (A17, A18) in Ref. pi) 



S Q (q , q; s) = V Q (q , q) + J — V Q (q , k) ^_ j p - ^ - ^— — S Q (k, « -) 

(14) 



Here the momentum Q represents the momentum of the pair with respect to the background, 
and Q the Pauli blocking operator which projects onto unoccupied intermediate states. For 
the above expression angle averages have been used, i.e. ||Q + k| 2 « iQ 2 + k 2 , and the 
s invariant has been defined as s = AE\^ + ^ — Q 2 . This approach, in the context of the 
Dirac-Brueckner-Hartree-Fock approximation, has been applied with reasonable success to 



the study of infinite nuclear matter [28] as well as finite nuclei ground state properties [29 



An appealing feature of the above equation for Q is its direct connection with the bare 
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NN potential model in free space. Indeed, adopting the same definitions as in Eqs. ([11]) and 
(0), 



, M . / Af . . 

£q(q,q; y/s) = ^/— £ Q (q,q;s)W — , (15) 



and 



V ti>4 = \I^Vtf,q)j¥-, (16) 

the following equation for the g-matrix is obtained 

( < r\ t/ / / \ _i_ f f f ff k \ MQ(Q,k) 

PQ(q , q V^) = V Q (q , q) + y ^V Q (q , k) ^— j i g _[_ |Q2 _ E| Q+k + fe * ^ 

(17) 

The above Dirac-Brueckner approach differs in a non-trivial way from the conventional 
nonrelativistic Brueckner approach. The density dependence of the one-boson-exchange 
interaction by means of effective Dirac spinors and the explicit relativistic kinematics are 
features which have no counterpart in the traditional Brueckner approach. However, with the 
suppression of these relativistic dynamical effects we can extract a minimum of relativistic 
features needed in the nonrelativistic model. In our approach selfconsistency is demanded 
with the following choice for the quasi-particle spectrum, 

El = p 2 + (M + U(p)) 2 . (18) 

As in the usual Brueckner-Bethe-Goldstone approach, the quasi-potential U(p) is obtained 
selfconsistently with the use of the continuous choice at the Fermi surface. As a first check, 
we consider the case of two nucleons interacting in free space, for which the Pauli blocking 
operator becomes the identity and the nuclear selfconsistent field vanishes. Furthermore, 
if the interaction is described in the pair cm. (Q = 0) then the (7-matrix corresponds to 
the free scattering matrix T as described by the BbS equation (cf. Eq. [13|). This limit 
is immediately verified upon substitution of y/s by 2\J q% + M 2 , with q a the on-shell cm. 
relative momentum. Thus, 

/H 3 k M 
(2^ (q ' k) go2 _ fc2 + ? / (k, q; Js) , (19) 
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which is the nonrelativistic Lippmann-Schwinger equation with the pole at the relativistically 
correct momentum. 

The other case of interest is infinite nuclear matter, where the relativistic structure of 
the quasi-particle spectrum introduced in Eq. QT8P can be assessed considering the following 
power expansion in terms of U/M 

E 2 JM « 2 + U(p)j + M(l + 0[{U/Mf}) . (20) 

If we substitute p 2 by the angle averaged quantity ^Q 2 + k 2 and proceed similarly with s, 
the similarity of the energy denominator in Eq. (O) with the one obtained using the nonrel- 
ativistic propagator is evident. An estimate of the accuracy of the above approximation at 
normal densities can be made considering (U/M) < 1/10. In such a case the above form of 
E 2 yields a propagator equivalent to its nonrelativistic counterpart with an accuracy better 
than 1%. This result supports the use of the nonrelativistic selfconsistent scheme to obtain 
the quasi-potentials. 

With the above considerations, the calculation of the two-body effective interaction 
needed in the full-folding optical potential proceeds with the proper choice of the s-invariant 
in Eq. fllTD , and consistent with the starting energy Q a defined in Eq. ([[]). Since Q a repre- 
sents the total pair energy of interacting nucleons with total momentum Q, the s-invariant 
is simply given by s = VL 2 a — Q 2 . 



V. APPLICATIONS 



A few remarks about the actual implementation of the computational procedures, aimed 
to obtain the full-folding optical potential U (E) , are worthy to note. First of all, we use the 
Slater approximation for the nuclear mixed density (f)'(j) in Eq. (JJ). This approximation has 
been discussed in the past || and becomes particularly suitable for the present applications, 
where we rely on the same point nuclear densities for the target ground state. Furthermore, 
we make use of the infinite nuclear matter ^-matrix to represent the effective interaction 
between the projectile and the target nucleon 0. As a result, the g-matrix evaluated at a 
Fermi momentum kp is folded with the target ground state density p at a local momentum 
k(R). The Slater approximation suggests the ansatz kp = k, where the local momentum k 
is determined from the density by k 3 (R) = 3ii 2 p(R)/2. These considerations yield for the 
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full-folding optical model potential 

U(k', k; E)=4irJ dR e l ( k '- k > R ( Pp (R) ^v(k', k) + p n (R) g nN (^', k)) , (21) 

where p p and p n are the local proton and neutron point densities respectively, and g n n rep- 
resent the off-shell Fermi-averaged amplitudes in the NN channel. For a particular channel 
this amplitude depends on the nuclear matter density via kp (implicit in g), and the local 
momentum k which sets the bounds for the Fermi motion of the target nucleons. This is 
expressed as 

g NN (k', k) = -4- / B(k - |P|)0K +P (k P ', k r ; Vi) dP , (22) 
Auk? J 

where K = (k + k')/2 and the relative momentum k r and k r ' are obtained following each 
of the two relativistic prescriptions discussed in Appendix A. We stress at this point that 
these amplitudes are calculated fully off-shell and that no assumption is made regarding 
the coordinate space structure of the <?-matrix. With these considerations the full-folding 
optical potential becomes a genuine non-local operator. Its use in the Schrodinger equation 
involves integro-differential equations which are solved exactly within numerical accuracy. 

The nuclear matter calculations for the g-matrix were done with fully self-consistent 
fields at various values of kp. We have considered the inversion potentials based upon 
the SPOO phase-shifts solution with Nijmegen-I, -II and Reid-93 reference potentials (spOO- 
NUl, SP00-NU2 and SP00-RE93 respectively), Argonne reference potential (sp00-av18), 
Paris reference potential (spOO-papj) and Gel'fand-Levitan-Marchenko inversion reference 
potential (spOO-invs). 

A. Medium and Fermi motion effects 

To disclose some of the features exhibited by the Fermi-averaged effective interaction in 
the context of the full-folding approach, we analyze the g averages in the proton-proton (pp) 
and neutron-proton (np) channels. Here we focus on the on-shell forward matrix element. 
In this case the g^N amplitude depends on the NA projectile momentum |k|, the Fermi 
momentum kp, and the local momentum k. In the case of the free t-matrix approach for 
the interaction we set Uf = but allow for the variation of P implied by the local nuclear 
density to account for the Fermi motion in the nucleus (|P| < k). In the case of a ^-matrix 
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element we set k = kp. In Fig. |2| we show the real and imaginary components of g pp and 
g np , based on the SP00-AV18 NNOMP, as functions of the beam momentum and for the 
sequence k = 0.6(0.2)1.4 fm _1 . In each frame we draw the case k = 0.6 fm -1 with a thick 
solid curve; the following values of k depart sequentially from this reference curve. The 
upper and lower frames correspond to g- and t-matrix results respectively. To assess the 
role of the imaginary part of the NNOMP, the results based on the full model are shown 
in the four leftmost frames with solid curves while the results where the imaginary part is 
suppressed are shown in the four rightmost frames with dashed curves. The two relativistic 
kinematics prescriptions -Giebink and AAY- yield almost indistinguishable results for the 
amplitudes. 

A comparison of results in the upper and lower frames indicates more dispersion due to 
the Fermi motion within the (/-matrix than in the t-matrix approach. This is a feature which 
is more pronounced in the real than the absorptive component of the amplitude. This is 
a clear indication of the role of the self-consistent fields in the Fermi-averaged quantities. 
In the context of the full NNOMP however, the manifestation of this sensitivity becomes 
diminished at projectile energies above 500 MeV as the real component of all amplitudes 
for |k| above ~5 fm _1 are significantly smaller than their imaginary counterparts. Such is 
not the case at the lower energies, where both the real and imaginary components of the 
amplitude become comparable. Thus, we do expect more sensitivity under Fermi motion in 
the context of the g matrix at projectile energies below 500 MeV. 

Another feature that emerges from Fig. |2] is the sensitivity of the absorptive part of g 
to the presence of the imaginary part of the NNOMP. Indeed, when this part is suppressed, 
the absorptive component of the amplitude saturates above ~5 fm -1 , in contrast with the 
full NNOMP where the trend of this absorption is to increase. The manifestation of this 
feature becomes clear when the Fermi average enters fully off-shell in the evaluation of the 
full-folding potential, as discussed in the following section. 

B. Total cross sections 

A global assessment of the full-folding model to the inclusion of relativistic kinematics 
and features of the underlying NN potential model over a wide energy range is made first 
by studying total cross sections for neutron-nucleus elastic scattering. In Fig. [3] we show 
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the measured |30] and calculated total cross sections for neutron elastic scattering from 16 0, 
40 Ca, 90 Zr and 208 Pb at beam energies ranging from 100 MeV up to 1 GeV. These cross 
sections are obtained by solving the scattering equations using full-folding optical potentials 
calculated as in Ref. ||, modified to include the kinematics discussed in this work. For 
completeness in this comparison we have considered the full NNOMP within the g- (solid 
curves) and i-matrix (dashed curves) approaches. In order to assess the role of the absorptive 
contribution of the NNOMP, we have also included results suppressing the imaginary part. 
These results are shown unmarked, whereas those using the full NNOMP have been labeled 
with a triangle at their right end. The data are represented with open circles. From this 
figure we observe a remarkable agreement between the full NNOMP ^-matrix results and 
the data, particularly at energies above ~200 MeV. Such is not the case for the t-matrix 
approach, or when the imaginary part of the bare NN potential is suppressed. In the former 
case the lack of nuclear medium effects becomes pronounced for projectile energies as high 
as 500 MeV in the case of Pb, but less for lighter nuclei. Conversely, the imaginary part of 
the NNOMP is crucial for the adequate description of cross section data at beam energies 
above 400 MeV, as observed when comparing the labeled and unlabeled curves. 

Total reaction cross sections for proton-nucleus elastic scattering at high energies are 
of increasing interest, particularly with current trends using spallation facilities with high 
energy beams. Thus, we calculated total reaction cross-section for elastic proton scattering, 
at beam energies up to 1.5 GeV, and compare them with available data pif . In strict 
analogy with the previous applications, and use of notations as in Fig. ||, in Fig. | we show 
g- (solid curves) and t-matrix (dashed curves) results wich used the full or suppressed the 
imaginary part of the NNOMP. Here again we have marked with triangular labels those 
results based on the full NNOMP, and left unmarked the ones with the imaginary part 
suppressed. The data are represented by open circles. The results shown in Fig. f| used 
the SP00-AV18 NNOMP. For this particular observable and energies of our study, we find 
very similar results among all the NNOMP i.e. irrespectively of the reference potential. 
Again we find little differences between the g- and t-matrix approaches above ~500 MeV. 
Thus, medium effects in the interaction are rather weak (albeit not neglegible) at these 
higher energies. Such is not the case below 400 MeV, particularly in the case of Pb, where 
a clear departure of the t-matrix with respect to the ^-matrix results is observed. However, 
these differences are smaller than the ones due to the presence of the imaginary part of the 



16 



NNOMP. In fact, we can see clear differences within the g- matrix approach by including and 
suppressing the imaginary part of the NNOMP at energies above 500 MeV. Since reaction 
cross section data are scarce above 700 MeV (only two data points for all four targets), the 
curves (marked with triangles) constitute a high energy prediction of our work. 

Another feature observed from Figs. ^ and |] is the trend of all cross sections to reach a 
plateau above 650 MeV. We have scrutinized more closely this feature and find that both 
o"t and or exhibit an almost linear dependence with A 2 /' 3 , with target masses 16< A <208. 
To illustrate this point we show in Fig. [5] the calculated ot and or as function of A 2 / 3 . The 
solid curves represent results from ^-matrix full-folding results, based on the SP00-AV18 full 
NNOMP, at various energies; the dotted curves represent the corresponding results without 
the imaginary part of the NNOMP. The dashed curves serve as reference and correspond to 
the parametric forms or = (—0.19 + O.lOA 2 / 3 ) and a T = (0.03 + 0.052v4 2//3 ) in barn units. 
Notice, only the full NNOMP results feature a moderate deviation from the parametric forms 
shown, which is not the case in calculations where the imaginary part of the NNOMP was 
suppressed (dotted curves). Thus, the absorptive part of the NNOMP inhibits the energy 
dependence of the total cross sections above 700 MeV. 

C. Differential cross sections 

Differential cross sections and spin observables for elastic scattering remain a challenge 
for any microscopic theory. An assessment of alternative relativistic kinematics is made 
with a comparison of these observables using Giebink and AAY relativistic approaches. 
These results are illustrated in Fig. |], where the differential cross section da/dQ, analyzing 
powers A y and spin rotation functions Q for p+ 40 Ca elastic scattering at 1.04 GeV are 
shown as function of momentum transfer q. We have used ^-matrices together with Giebink 
(solid curves) and AAY (dashed curves) kinematics. This figure shows the case where the 
differences should be most pronounced. Indeed, all the other cases show nearly complete 
overlap between the two cases. Nevertheless, Fig. [(] shows that the two prescriptions yield 
similar results for all observables. Some differences can be seen in A y , at momentum transfers 
above 1.5 fm _1 , where the AAY lie slightly above the Giebink kinematics results. We selected 
for all further calculations the Giebink kinematics. 

Next, we study six NNOMP in NA elastic scattering at various energies with the full- 
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folding model using the g-matrix approach. We have chosen scattering of protons from 40 Ca 
and 208 Pb for which there is a large body of high precision data over a wide energy range 
f32fl . The applications used the full NNOMP with reference to SPOO-NUl, -nij2, -pari, 
-AV18, -RE93 and -invs solutions. As in most cases the differences among these reference 
NN potentials are moderate, we have chosen all curves with a single pattern. This also helps 
to illustrate the level of sensitivity of NA scattering upon on-shell equivalent potentials with 
different off-shell behavior. 

Considering that the results presented here correspond to parameter free calculations, the 
overall description of all elastic scattering data is remarkably good. Only a close examination 
of the results shows the limitations of our approach. In Fig. |7] for p+ 40 Ca scattering at 
300, 400 and 497.5 MeV, for instance, we observe a very good agreement between the full- 
folding model results for all observables and momentum transfers above ~1 fm -1 . However, 
some discrepancies remain in the description of low momentum transfer data. In Fig. |8| 
are shown p+ 40 Ca results for energies between 650 and 1040 MeV. We observe two curves 
slightly separated from the rest. Such curves correspond to the SP00-AV18 and spOO-pari 
NNOMP, affecting mainly the differential cross section and analyzing powers at 800 and 
1040 MeV. The possible cause of this feature, and the fact that this is more pronounced in 
the case of Ca than in Pb is not fully understood but numerical reasons cannot be ruled out. 
A similar trend is observed in Fig. ^ for p+ 208 Pb scattering. These limitations have already 
been noticed in previous nonrelativistic full-folding calculations. Thus, the identification of 
a unique missing element in the approaches, to account for the low q failures remains an 
open issue. Some possibilties have been discussed elsewhere ||. 

In Figs. H and |10| we show results for p+ 208 Pb scattering at energies between 300 MeV 
and 1 GeV. With the exception of the 400 and 497.5 MeV cases, the cross section is very 
well described by our calculations. With respect to the spin observables, there is a tendency 
to lose structure relative to the measured values at energies above 650 MeV. 

D. Approximations 

The calculation of optical potentials, within the full-folding approach used herein, are 
made without refuge in assumptions either about the local structure of the NN effective 
interaction or about the final structure of the NA coupling. In fact, these potentials are 
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treated as non-local operators and are the result of a detailed account of the NN effective 
interaction off-shell. In contrast, the optical potentials in the nonrelativistic impulse approx- 
imation, discussed in Ref. ||, are assumed local and constructed with only on-shell t-matrix 
elements as effective interaction. Important differences have been observed between these 
two approaches when applied to intermediate energy (200-400 MeV) NA scattering. An 
illustration of these differences is made in Fig. |TT| for 800 MeV p+ 208 Pb elastic scattering, 
where the differential cross section, analyzing powers and spin rotations are shown as func- 
tion of the momentum transfer. The g-matrix full-folding results are represented with solid 
curves. The t-matrix full-folding results are represented with long-dashed curves, whereas 
short-dashed and dotted curves are used to represent the results of the off-shell tp and 
on-shell tp results respectively. Clear differences can be seen between the g- and t-matrix 
results, particularly for the spin observables. These are pronounced for q > 1.5 fm _1 which 
illustrates the level of sensitivity to medium effects at these momentum transfers and for 
this particular case. Within the same t-matrix approach however, the full-folding results are 
similar to those obtained within the tp approximation. The extent of these sensitivities is 
comparable to contributions from short range correlations ||. 

The full-folding calculations presented here are first of their kind to be tested at energies as 
high as ~1.5 GeV. Differences are observed between these and previous results, particularly 
in the 800 MeV applications at forward angles H. Apart from the locality issue in the NA 
coupling and the treatment of the NN effective interaction off-shell, additional effects have 
been included in those local potential calculations and which may cause variations between 
the results. Particularly relevant seem to be the short range correlations and electromagnetic 
spin-orbit contributions to the NA coupling. A careful assessment of these effects in the 
context of the full-folding approach is needed. Considering that calculations in the tp scheme 
are much simpler than those in the full-folding approach, the former becomes quite suitable 
for exploratory purposes. Quantitative comparisons, however, do require the inclusion of 
medium effects within the full-folding approach. 

VI. SUMMARY AND CONCLUSIONS 

The nonrelativistic full-folding optical model approach for NA scattering has been ex- 
tented into the relativistic regime. Kinematical issues involving the off-shell Lorentz boost 
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of the colliding particles between the NN and the NA cm. reference frames have been 
addressed. Also, explicit nuclear medium effects have been incorporated with the use of mi- 
croscopic NN effective interaction as obtained in the framework of nuclear matter (/-matrix 
using an NNOMP which fully accounts for the inelasticities and isobar resonances at nu- 
cleon energies as high as ~3 GeV. The nuclear matter (/-matrices were obtained considering 
both Pauli blocking and selfconsistent nuclear fields as in the traditional Brueckner theory. 
Minimal relativity corrections were extracted from the Brockmann and Machleidt approach 
to relativistic nuclear matter. Effects arising from Wigner rotations and electromagnetic 
spin-orbit corrections were not included. 

The study considered both t-matrix and in-medium self-consistent ^-matrix approaches. 
With the inclusion of relativistic kinematics corrections, in conjunction with a realistic de- 
scription of NN resonances and inelasticities by means of an NNOMP, we obtain a good 
description of both the total and the differential scattering observables for NA elastic scat- 
tering. The results exhibit a weak sensitivity to the choice of the relativistic approach -AAY 
or Giebink- to correct the kinematics. We also observe that medium effects are significant 
over the whole energy range of our study, although they are rather weak above 400 MeV. In 
contrast, the inelasticities of the NN interaction become important above 400 MeV as was 
observed particularly in the description of the total cross sections. 

Although our study allows a reasonable description of the differential observables at 
energies as high as 1 GeV in NA scattering, specific details remain to be explained; notably 
the low q behavior of the spin observables in the 400-500 MeV range. As our primary 
effort has been to provide a parameter free nonrelativistic framework for the study of NA 
elastic scattering with a minimal account of relativistic effects at these high energies, a 
systematic study of various other effects has not been pursued in order to keep our discussion 
focused. However, future work will scrutinize the systematics of the calculated observables 
under the use of alternative densities, mixed density representations, electromagnetic effects 
and higher order correlations. Off-shell effects arising from the non-locality range of the 
separable description of the NNOMP above pion threshold also require further investigation. 
Nevertheless, the level of agreement we have achieved, within the nonrelativistic approach, 
is comparable to what has been obtained within relativistic approaches. 

The study presented is limited in its value as we have not used a covariant two-body and 
(A+l)-body dynamics. The approach adopted is a practical one and is largely motivated 
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and justified by the good agreement of the numerical results with data. Thus we claim, 
that use of minimal relativity in conjuction with NN optical model interactions, which 
fully account for the inelasticities and isobar resonances above the pion threshold, yield 
quantiative descriptions of NA scattering up to energies of 1.5 GeV. 



APPENDIX A: RELATIVISTIC KINEMATICS 



We have used two prescriptions for the relativistic kinematics involved in the transfor- 
mation of the two-body colliding momenta between the pair cm. and the projectile-nucleus 
cm. These schemes have been developed by Giebink |23[ and by Aaron et al. pif] . The lat- 



ter has also been obtained by imposing time reversal invariance on the scattering amplitude 



Following Giebink in the context of a manifestly Lorentz invariant two-body transition 
amplitude, let k = (a), k) andp = (e, p) be the projectile and struck-nucleon four-momenta in 
the projectile-nucleus cm. The corresponding momenta in the exit channel are represented 
with k! and p' . In Giebink's approach the total four-momentum is conserved, 

k + p = k' + p' = Q, (Al) 

which, to be valid, requires each particle to be off mass-shell. The relative momenta in the 
two-body cm. is readily obtained by applying a Lorentz transformation to momenta with 
a single boost (3. This boost is obtained from the invariant 

(cu + e) 2 -(k + p) 2 = S , (A2) 

and from 

k_fp_k^+jy 

With this velocity, a direct calculation yields for the incident (k r ) and outgoing (k' r ) relative 
momenta 

k _ (e + e r )k-(u + u r )p k ,_ (t + e' r )k'-(u)' + u' r )p' 
e + e r + uj + u r ' r ~~ e' + e' r + Q' + w' r 

In these expressions the subscript V denotes the on mass-shell relative energy 



u r = \l m 2 p + k r , e r = \jm\ + k r , (A5) 
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where m p and m t represent the masses of the projectile and struck-nucleon respectively. In 
all the above expressions the magnitude of the relative momentum k r 2 is needed. It can be 
shown that 

k r 2 = ^m^P 2 ), (A6) 

where the ^-function is defined as 



y, z) = y/(x - y - zf - Ayz . (A7) 

Notice that the ^-function in Eq. (|A6|) is evaluated at the off- mass-shell invariants k 2 = uj 2 — 
k 2 and p 2 = e 2 — p 2 . The actual implementation of Giebink's procedure faces the difficulty 
of improper Lorentz transformations which occur when S in Eq. (|A2|) becomes negative; 
this happens for very large momenta k + p. However, as the bound-state wavefunctions of 
the struck nucleons confine the momentum distribution of p to magnitudes below ~2 fm _1 , 
such improper contributions occur for very large k and k', affecting only far off-shell matrix 
elements in C/(k',k). A way to circumvent this difficulty is to restrict S in Eq. ( |A2j) near 
on-mass-shell. Thus, we approximate 

S » m\ + m 2 + 2ue - 2k ■ p . (A8) 

Averaging the Fermi motion of the target nucleons (p) allows the simplification 

S -> S = m 2 + m\ + 2cue . (A9) 

Consistent with the above, the following forms are obtained for the relative energies 

u)E + m 2 , ooe + m 2 
uj r = and e r = -=-^ . (A10) 

V 5*0 V S 



Clearly vSo = uj r + e r , which is a handy result for the energy denominators in Eq. flA4| 
The full-folding calculations following Giebink relativistic kinematics were made using e ■■ 
e' = M and Q — u' — M + E^i,, with M the nucleon mass. 



In the relativistic prescription developed by Aaron and collaborators [^TJ, the relative 
momenta for the incoming and outgoing channels require different boost velocities. In each 
channel the particles are set on-mass-shell and the corresponding boost is represented by 



Eq. (|A3|) . The resulting relative momentum exhibits the same structure as the one given 



22 



by Eq. fl^j|) but with the substitution uj — > c<j(k) and £ — > e(p). A direct calculation yields 
for the magnitude 

kl = ^-e(s tn ,m 2 p ,m 2 t ) , (All) 

where 

Sm = (e(p)+o;(k)) 2 -(p + k) 2 . (A12) 
An analogous result is obtained for the outgoing channel. 

APPENDIX B: ALGORITHM TO DETERMINE THE NNOMP FROM DATA 

Consider that there are three distinct Hamiltonians They are the reference Hamil- 
tonian H , a projected Hamiltonian Hpp, and a full optical model Hamiltonian TC. The 
first of these, the reference Hamiltonian Hq := T + Vq, invokes a given potential Vo for 
which one can find Schrodinger equation reference solutions. The physical outgoing solu- 
tions ip := i/Jq(t, k, E) of H we suppose yields a unitary S- matrix. We assume further that 
this Hamiltonian is completely specified such that evaluation of any quantity, wave function, 
S-matrix, K-matrix etc. is facilitated. The Feshbach projection operator formalism gives 
the projected Hamiltonian, PH P = H PP . We presuppose completeness, P + Q — 1, and a 
finite rank representation of the Q space 

N N 

Q:=^2\9i><9i\ = ^2\i><i\, (Bl) 

i=i i=i 

with the Q space basis functions |<&$ > interpreted as doorway states. With these doorway 
states we make the link between the QCD and the hadronic sectors; the latter encompassing 
nucleons, mesons and other free particles. Thus we assume that meson creation/annihilation 
occurs only in the highly nonlinear QCD sector so that Q space wave functions are projections 
of such processes onto hadronic particle coordinates. The third of our Hamiltonians, the full 
optical model Hamiltonian, comprises the reference Hamiltonian H and the proper optical 
model potential V. That potential is complex and nonlocal, and taken to be separable of 
finite rank, H := T + Vo + V(r, r'\ Isj, E). 

The Schrodinger equation specified with 7i has regular physical solutions : = 
\l/ + (r, k, E) whose asymptotic boundary conditions we deem to match with the experimental 
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elastic channel S-matrix. Specifically, for these experimental S-matrices we have used the 
continuous solutions SPOO from VPI/GWU |27[]. The reference potential Vo and separable 
potential form factors are still to be specified in detail with the application. 

To obtain the optical potential on the basis of a given reference potential, we express first 
the solutions of the projected Hamiltonian in terms of the reference Hamiltonian and the a 
priori defined Q space projector. The Lippmann-Schwinger equation, 

|Vp >= 1-00 > ~J2 G o\j >< J\ h Qp\^p >> (B2) 
5 

is still very general and does not depend upon a specific representation. However, in the 
following we assume a partial wave expansion and the following equations are identified as 
radial equations with the set of quantum numbers suppressed. 
As projector orthogonality PQ = QP = implies that 

N 

< j\H QP \4> P >= J2i< *\Gt\$ >}j- < #o >, (B3) 

% 

the solutions of Eq. flB^) can be written in terms of |^ > as 

N N 

|Vp > = |^o >-J2 G o\i > {< $ l G o 1$ >}? < Mo > = lV>o > -J2 G o A iMo >, 

ij ij 

(B4) 

wherein one can identify a separable potential 

K>{<$|G + |$>} i r. 1 <j| = \i>\ ij <j\=:A ij (r,r'). (B5) 

Note then that definition of Q space gives a specification of the separable strengths Xij(lsj, E) 
that is unique. The resultant Eq. ( |B4| ) has the form of a first order Born approximation but 
in fact it is an exact result. 

To proceed, initially we abandon the exactitude of Eq. (fB3j) and require the strength 
matrix, 

Ay = {< $1(7+1$ >}r.\ (B6) 

to be constrained asymptotically by the experimental S-matrix of the full Hamiltonian 
Schrodinger equation, i.e. asymptotically we induce \ipp >= >■ This implies that 
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complex optical model strengths A^ emerge as a result of matching to Riccati-Hankel func- 
tions and non unitary S-matrices with 

\y n >= |^p >^ (-h~(rk) + h + (rk)S(k)) . (B7) 

The strengths \j then can be simply determined from the linear system of equations 

^h + (Rk) (S(k) - S (k)) = J^G+\i > *U < Mo >■ (B8) 

ij 

To reinforce a Lippmann-Schwinger equation, with the experimental S-matrix as boundary 
condition or equivalently with strengths Ay from Eq. (|B8|), a transformation of the separable 
potential Eq. (|B5| ) is made. This is achieved with 

v(ry): = A (i^W (B9) 

which contains the separable potentials as defined with Eq. (p5|) but whose strengths now are 
solutions of Eq. (|B8|). As the transformation Eq. ( |B9| ) contains integration of orthonormal 
functions, only strengths are altered. Using this optical model in the full Hamiltonian, 
physical solutions are obtained with reference solutions \ipo > and Greens function Gq of 
the reference Hamiltonian Hq by means of the Lippmann-Schwinger equation 

\*h >= l^o > +G$V\*h > . (BIO) 
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FIG. 1: Schematic representation of the collision of the projectile with a target nucleon. Quantities 
in parentheses represent the four-momenta of the colliding nucleons. 
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FIG. 2: The Fermi-averaged forward amplitude (jnn at &f = 1 fm _1 (upper frames) and tjvTV 
(lower frames) based on the sp00-av18 NNOMP at local momenta £=0.6(0.2)1.4 fm _1 in the 
pp and np channels, and as functions of the projectile momentum. The four leftmost frames 
represent results using the full NNOMP whereas the four rightmost frames with the dashed curves 
correspond to results whith the imaginary part of the NNOMP suppressed. The thick solid curves 
denote results for k=0.6 fm _1 . 
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FIG. 3: Total cross section for neutron elastic scattering from Pb, Zr, Ca and O as 



functions of the projectile energy. The data [3C] are represented with open circles. The solid 
and dashed curves represent full-folding results using the g- and t-matrix respectively. The curves 
corresponding to the full NNOMP are marked with a triangular label at their right end, whereas 
those results with the imaginary part of the NNOMP suppressed are unmarked. 



31 



~ i — i — i — r 



U 



o 



] 



_1 I I u 



_l \/\ I I /I / [_ 



o 
o 



o 
o 

CM 



o 
o 

O 



o 
o 

00 



o 
o 



o 
o 



- o 

— o 

- cs 



in 



d 



o 
d 



[q] *o 



> 



>1 



FIG. 4: Predicted reaction total cross section for proton elastic scattering from 208 Pb, 90 Zr, 40 Ca 
and 16 as functions of the beam energy. The data were taken from Ref. [31]. The solid and 
dashed curves represent full-folding results using the g- and £-matrix respectively. The curves 
corresponding to results based on the full NNOMP have been marked with a triangular label at 
their right end, whereas those results with the imaginary part of the NNOMP suppressed are 
unmarked. 



32 




FIG. 5: The reaction cross section (<jr) and total cross section (<tt) for proton and neutron elastic 
scattering as function of A 2 ! 3 for projectile energies of 650, 800, 1040, 1250 and 1500 MeV. The 
solid and dotted curves correspond to ^-matrix full-folding results based on the full and imaginary- 
suppressed NNOMP respectively. The dashed curves correspond to straight lines in terms of A 2 ! 3 
(see text for details). 
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FIG. 6: Sensitivity to the kinematics: the calculated differential cross-section (upper frame), 
analyzing power (middle frame) and spin rotation (lower frame) as functions of the momentum 
transfer for p+ 40 Ca elastic scattering at 1.04 GeV. The results correspond to nonrelativistic in- 
medium full-folding optical potentials using Giebink's (solid curves) and AAY (dashed curves) 
relativistic kinematics. 
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FIG. 7: Calculated differential cross-section (upper frames), analyzing power (middle frames) and 
spin rotation (lower frames) as functions of the momentum transfer for p+ Ca elastic scattering 
at 300 MeV, 400 MeV and 497.5 MeV. All curves are obtained from in-medium full-folding optical 
potentials using relativistic kinematics. The data are from Ref. [p 7 ^]. 
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FIG. 9: Calculated differential cross-section (upper frame), analyzing power (middle frame) and 
spin rotation (lower frame) as functions of the momentum transfer for p+ 208 Pb elastic scattering 
at 300 MeV, 400 MeV and 497.5 MeV. All curves are obtained from in-medium full-folding optical 
potentials using relativistic kinematics. The data are from Ref. [p 7 ^]. 
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FIG. 11: A comparison of (/-matrix full- folding (solid curves), free t-matrix full- folding (long- 
dashed curves) off-shell 't/?' (short-dashed curves) and on-shell 'tp' (dotted curves) results in the 
case of p+ 208 Pb elastic scattering at 800 MeV. 
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